Using the Scaling Analysis to Characterize 
Financial Markets 



T. Di Matteo a ' b , T. Aste b , M. M. Dacorogna '* 

a INFM - Dipartimento di Fisica "E. R. Caianiello" , Universitd degli Studi di 
Salerno, 84081 Baronissi (SA), Italy. 

b Applied Mathematics, Research School of Physical Sciences, Australian National 
University, 0200 Canberra, Australia. 

c Converium Ltd, General Guisan - Quai 26, 8022 Zurich, Switzerland. 



Abstract 

We empirically analyze the scaling properties of daily Foreign Exchange rates, Stock 
Market indices and Bond futures across different financial markets. We study the 
scaling behaviour of the time series by using a generalized Hurst exponent approach. 
We verify the robustness of this approach and we compare the results with the 
scaling properties in the frequency-domain. We find evidence of deviations from the 
pure Brownian motion behavior. We show that these deviations are associated with 
characteristics of the specific markets and they can be, therefore, used to distinguish 
the different degrees of development of the markets. 
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1 Introduction 



The scaling concept is increasingly applied outside the traditional physical 
sciences domain [1,2,3,4,5,6,7]. In the recent years, its application to financial 
markets, initiated by Mandelbrot in the 1960 [8,9], has largely increased also 
in consequence of the abundance of available data [1]. Two types of scaling 
behaviors are studied in the finance literature: 

(1) The behavior of some forms of volatility measure (variance of returns, 
absolute value of returns) as a function of the time interval on which the 
returns are measured. (This study will lead to the estimation of a scaling 
exponent related to the Hurst exponent.) 

(2) The behavior of the tails of the distribution of returns as a function of the 
size of the movement but keeping the time interval of the returns con- 
stant. (This will lead to the estimation of the tail index of the distribution 
[7]-) 

Although related, these two analysis lead to different quantities and should not 
be confused as it is often the case in the literature as can be seen in the papers 
and debate published in the November 2001 issue of Quantitative Finance. For 
more explanations about this and the relation between the two quantities, the 
reader is referred to the excellent paper by Groenendijk et al. [10]. In this study, 
we are interested in the first type of analysis. Until now, most of the work has 
concentrated in studies of particular markets: Foreign Exchange [1,7,11], Stock 
[12] or Fixed Income [13]. These studies showed that empirical scaling laws 
hold in all these markets and for a large range of frequencies: from few minutes 
to few months. In a recent book [7] , the hypothesis of heterogeneous market 
agents was developed and backed by empirical evidences. In this view, the 
agents are essentially distinguished by the frequency at which they operate 
in the market. The scaling analysis, which looks at the volatility of returns 
measured at different time intervals, is a parsimonious way of assessing the 
relative impact of these heterogeneous agents on price movements. Viewing the 
market efficiency as the result of the interaction of these agents [14], brings 
naturally to think that it is the presence of many different agents that would 
characterize a mature market, while the absence of some type of agents should 
be a feature of less developed markets. Such a fact should then reflect in the 
measured scaling exponents. The study of the scaling behaviors must therefore 
be an ideal candidate to characterize markets. To further explore this issue, we 
perform an empirical analysis of daily data across different financial markets 
to examine the similarities or differences in the scaling properties. 

Until the 1960s, the only stochastic and scaling model in finance was the 
Brownian motion, originally proposed by Bachelier in 1900 [6,15], and de- 
veloped several decades later [16]. This theory predicts that the returns of 



2 



market prices should follow a normal distribution with stable mean and finite 
variance. However, there are ample empirical evidences that the returns are 
not normally distributed but have higher peak around the mean and fatter 
tails [6,7]. Moreover, it is also observed that volatility clustering is a general 
characteristic of financial markets [7]. Generalizations of the classical Brow- 
nian motion were made by Mandelbrot and followers involving either frac- 
tional Brownian motions [9,17,18], or Levy motion [8,19,20,21,22,23]. Closely 
related additive scaling models have also been developed [6,12,24]: Brownian, 
fractional Brownian, Levy processes. The above approaches generally involve 
additive monofractal processes and analyses; but, in contrast, several scal- 
ing systems appear to be more complex. Recently, a controversy has erupted 
[25,26,27,28] to know if the processes that describe financial data are truly 
scaling or simply an artifact of the data. Moreover, several publications pro- 
pose new scaling models or empirical analyses that better describe empirical 
evidences [25,26,27,28,29]. It should be however noted that - as underlined by 
Stanley et al. [30] - in statistical physics, when a large number of microscopic 
elements interact without characteristic scale, universal macroscopic scaling 
laws may be obtained independently of the microscopic details. 

In this paper we address the question of the scaling properties of financial time 
series by empirically analyzing daily data for Foreign Exchange rates, Stock 
Market indices and Bond futures (described in Section 2). We study very de- 
veloped as well as emerging markets in order to see if the scaling properties 
differ between the two and if they can serve to characterize and measure the 
development of the market. Here the scaling law is not used to conclude any- 
thing on the theoretical process but on the contrary we use it as a "stylized 
fact" that any theoretical model should also reproduce. The purpose of this 
paper is to point out how a relatively simple statistics gives us indications on 
the market characteristics, very much along the lines of the review paper by 
Brock [31]. In Section 3, we recall the theoretical framework and we intro- 
duce the generalized Hurst exponents analysis. In Section 4 we describe the 
methodology utilized to empirically analyze the data. In Section 5 we compute 
and compare the scaling spectral exponents and the Hurst exponents. Finally 
some conclusions are given in Section 6. 



2 The Studied Markets 



We study several financial markets which are at different development stage: 
mature and liquid markets, emerging and less liquid markets. These markets 
deal with different instruments: equities, foreign exchange rates, fixed income 
futures. In particular, the data that we analyze are: Foreign Exchange rates 
(FX) (see Table 1), Stock Market indices (SM) (see Table 2), Treasury rates 
corresponding to twelve different maturity dates (TR) (see Table 3) and Eu- 



3 



rodollar rates having maturity dates ranging from 3 months to 4 years (ER) 
(see Table 4). Hereafter we give a brief description of the time-series studied 
in this paper. 

• FX: The Foreign Exchange rates (Table 1) are daily rates of exchange of 
major currencies against the U.S. dollar. The time series that we study go 
from 1990 to 2001 and 1993 to 2001. These rates have been certified by 
the Federal Reserve Bank of New York for customs purposes. The data are 
noon buying rates in New York for cable transfers payable in the listed cur- 
rencies. These rates are also those required by the Securities and Exchange 
Commission (SEC) for the integrated disclosure system for foreign private 
issuers. The information is based on data collected by the Federal Reserve 
Bank of New York from a sample of market participants. 

• SM: The Stock Market indices (reported in Table 2) are 32 of the major 
indices of both very developed markets and emerging markets. These daily 
time series range from 1990 or 1993 to 2001. 

• TR: The Treasury rates (Table 3) are daily time series going from 1990 
to 2001. The yields on Treasury securities at 'constant maturity' are in- 
terpolated by the U.S. Treasury from the daily yield curve. This curve, 
which relates the yield on a security to its time to maturity, is based on the 
closing market bid yields on actively traded Treasury securities in the over- 
the-counter market. These market yields are calculated from composites of 
quotations obtained by the FD Bank of New York. The constant maturity 
yield values are read from the yield curve at fixed maturities, currently 3 
and 6 months and 1, 2, 3, 5, 7, 10, and 30 years. The Treasury bill rates 
are based on quotes at the official close of the U. S. Government securities 
market for each business day. They have maturities of 3 and 6 months and 
1 year. 

• ER: The Eurodollar interbank interest rates (Table 4) are bid rates with 
different maturity dates and they are daily data in the time period 1990- 
1996 [32]. 

As an example, the behaviors of FX rates and the SM rates for Japan ( JPY/USD 
and Nikkei 225) and Thailand (THB/USD and Bangkok SET) as a function 
of time t are shown in Fig. 1 in the time period 1997-2001. Another example 
is given in Fig. 2 that shows the TR time series as a function of t at different 
maturities dates in the time period 1997-2001 and the ER time series as a 
function of t at different maturities dates in the time period 1990-1996. 



3 The theoretical framework 

The scaling properties in time series have been studied in the literature by 
means of several techniques such as the rescaled range statistical analysis 
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R/S [33] and modified R/S analysis [34], multiaffine analysis, detrended fluc- 
tuation analysis (DFA) [35,36,37,38,39,40], periodogram regression (GPH method) 
[41], the ARFIMA estimation by exact maximum likelihood [42,43], the (m, k)- 
Zipf method [44], the moving-average analysis technique [45], the Average 
Wavelet Coefficient Method [46,47,48,49], and -in this respect- connection to 
multi-fractal/multi-affine analysis (the q order height-height correlation) have 
been made [50,51,52]. The rescaled range statistical analysis (R/S analysis) 
was first introduced by Harold Edwin Hurst [33] to describe the long-term 
dependence of water levels in rivers and reservoirs. It provides a sensitive 
method for revealing long-run correlations in random processes. This analy- 
sis can distinguish random time series from correlated time series and gives 
a measure of a signal "roughness". What mainly makes the Hurst analysis 
appealing is that all these information about a complex signal are contained 
in one parameter only: the Hurst exponent. The original approach of Hurst to 
the scaling properties of the time series is recalled in details in Appendix. One 
of the weaknesses of the original method is that it relies on maximum and 
minimum data, which makes it very sensitive to outliers. In order to study the 
multi-fractal features of the data we here use an alternative method to the 
original approach of Hurst. 



3. 1 Generalized Hurst exponent 



The Hurst analysis brings to light that some statistical properties of time 
series X(t) (with t=u, 2z/, ku, T) scale with the observation-period (T) 
and the time- resolution (u). Such a scaling is characterized by an exponent 
H which is commonly associated with the long-term statistical dependence 
of the signal. A generalization of the approach proposed by Hurst should 
therefore be associated with the scaling behavior of statistically significant 
properties of the signal. To this purpose we analyze the g-order moments of 
the distribution of the increments [9,53] which is a good characterization of 
the statistical evolution of a stochastic variable X(i). It is given by the: 



K (r) = (\ X (t + -)- X W q ) (1) 



where the time-interval r can vary between v and r max . (Note that, for q = 2, 
the K q (r) is proportional to the autocorrelation function a(r) = (X(t + r)X(t)).) 
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The generalized Hurst exponent H(q) 1 can be defined from the scaling be- 
havior of K q {r) which can be assumed to be given by the relation [53] 



Within this framework, we can distinguish between two kind of processes: (i) 
a process where H(q) = H, constant independent of q; (ii) a process with 
H(q) not constant. The first case is characteristic of uni-scaling or uni- fractal 
processes and its scaling behavior is determined from a unique constant H 
that coincides with the Hurst exponent H . This is for instance the case for 
self-affme processes where qH(q) is linear (H(q) = H) and fully determined 
by its index H. (Recall that, a transformation is called afline when it scales 
time and distance by different factors, while a behavior that reproduces itself 
under affine transformation is called self-affme [9] . A time-dependent self-affme 
function X(t) has fluctuations on different time scales that can be rescaled so 
that the original signal X(t) is statistically equivalent to its rescaled version 
\~ H X(\t) for any positive A [54], i.e. X(t) ~ X~ H 'X(Xt). Brownian motion 
is self- affine by nature.) In the second case, when H(q) depends on q, the 
process is commonly called multi-scaling (or multi-fractal) [54,55] and different 
exponents characterize the scaling of different g-moments of the distribution. 
In this 'curve' of exponents H(q), some values of q are associated with special 
features. For instance, when q — 1, H(l) describes the scaling behavior of the 
absolute values of the increments. The value of this exponent is expected to 
be closely related to the original Hurst exponent, H, that is indeed associated 
with the scaling of the absolute spread in the increments. The exponent at 
q = 2, is associated with the scaling of the autocorrelation function and is 
related to the power spectrum [56] . A special case is associated with the value 
of q = q* at which q*H(q*) = 1. At this value of q, the moment K q *(r) scales 
linearly in r [9]. Since qH(q) is in general a monotonic growing function of 
q, we have that all the moments H q (r) with q < q* will scale slower than 
r, whereas all the moments with q > q* will scale faster than r. The point 
q* is therefore a threshold value. In this paper we focalize the attention on 
the case q — 1 and 2. Clearly in the uni-fractal case H(l) = H(2) = H(q*). 
Their values will be equal to 1/2 for the Brownian motion and they would 
be equal to H ^ 0.5 for the fractional Brownian motion. However, for more 
complex signals, these coefficients do not in general coincide. We thus see that 
the non-linearity of the empirical function qH(q) is a solid argument against 
Brownian, fractional Brownian, Levy, and fractional Levy models, which are 
all additive models therefore giving for qH(q) straight lines or portions of 
straight lines. The curves for qH(q) vs. q are reported in Fig. 3 for some of 
the data. One can observe that, for all these time series, qH(q) is not linear 

1 We use H without parenthesis as the original Hurst exponent and H(q) as the 
generalized Hurst exponent. 




(2) 
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in q but slightly bending below the linear trend. The same behavior holds 
for the other data. This is a clear sign of deviation from Brownian, fractional 
Brownian, Levy, and fractional Levy models, as already seen in FX rates [1]. 
(Other cases showing marked deviations from Brownian motion have been 
discussed elsewhere [50,52,57,58,59].) 

3.2 Scaling spectral density and Hurst exponent 

For financial time series, as well as for many other stochastic processes, the 
spectral density S(f) is empirically found to scale with the frequency / as 
a power law: S(f) oc f' 13 . Here we use a simple argument to show how this 
scaling in the frequency domain is related to the scaling in the time-domain. 
Indeed, it is known that the spectrum S(f) of the signal X(t) can be conve- 
niently calculated from the Fourier transform of the autocorrelation function 
(Wiener-Khinchin theorem). On the other hand, the autocorrelation function 
of X(t) is proportional to the second moment of the distribution of the in- 
crements which, from Eq. 2, is supposed to scale as K 2 ~ t 2H ^ . But, the 
components of the Fourier transform of a function which behaves in the time- 
domain as r a are proportional to ]~ a ~ x in the frequency-domain. Therefore, 
we have that the power spectrum of a signal that scales as Eq. 2 must behave 
as: 

S(f) oc f- m *>- X ■ (3) 

Consequently, the slope (5 of the power spectrum is related with the general- 
ized Hurst exponent for q = 2 through: = 1 + 2H(2) [56]. Note that Eq. 3 
is obtained only assuming that the signal X(t) has a scaling behavior in ac- 
cordance to Eq. 2. Here we are not making any hypothesis on the kind of 
underlying mechanism that might lead to such a scaling behavior. 



4 Methodology 

Let us here recall that the theoretical framework we presented in the previous 
section is based on the assumption that the signal has the scaling property 
described in Eq. 2. Moreover, we have implicitly assumed that the scaling 
properties associated with a given time series stay unchanged across the ob- 
servation time window T. On the other hand, it is well known that financial 
time series show evidences of variation of their statistical properties with time, 
and depend on the observation time window T. The simplest case which shows 
such a dependence is the presence of a linear drift (r/t) added to a stochastic 



7 



signal (X(t) = X(t) +rjt) with X(t) satisfying Eq. 2 and the above mentioned 
properties of stability within the time window. Clearly, the scaling analysis 
described in the previous section must be applied to the stochastic component 
X(t) of the signal. This means that we must subtract the drift rjt from the 
signal X(t). To this end one can evaluate rj from the following relation: 

(X(t + r)-X(t)}=r ] r . (4) 



Other more complex deviations from the stationary behavior might be present 
in the financial data that we analyze. In this context, the subtraction of the 
linear drift can be viewed as a first approximation. The accuracy of this ap- 
proximation has been verified by varying the observation time window. We 
verify that the results obtained for different time window sizes are all compa- 
rable, with fluctuations within a range of 10%. 

Our empirical analysis is performed on the daily time series TR, ER, FX and 
SM (described in Section 2) which span typically over periods between 1000 
and 3000 days. In particular, we analyze the time series themselves for the TR 
and ER, whereas we compute the returns from the logarithmic price X(t) = 
ln(P(t)) for FX and SM. Moreover, all of these variables are 'detrended' by 
eliminating the linear drift as described above (Eq. 4). 

We compute the g-order moments K q {r) (defined in Eq. 1) of the 'detrended' 
signals and their logarithms with r in the range between v — 1 day and r max 
days. In order to test the robustness of our empirical approach, for each series 
we analyze the scaling properties varying r max between 5 and 19 days. The 
resulting exponents computed using these different r max are stable in their 
values within a range of 10%. 

We verify that the scaling behavior given in Eq. 2 is well followed and we 
compute the associated generalized Hurst exponent H(q) whose values are 
given in Section 5. 

In order to test that our method is not biased we estimate the generalized 
Hurst exponents for simulated random walks produced by using three differ- 
ent random numbers generators. We perform 100 simulations of random walks 
with the same number of data points as in our samples (991 and 3118) and 
estimate the generalized Hurst exponents H(l) and H(2) and the power spec- 
tra exponents (5. The results are reported in Table 5. In all the cases, H(l) 
and H{2) have values of 0.5 within the errors. Only when we consider uni- 
formly distributed random numbers in the interval (0,1) (Rand uses a lagged 
Fibonacci generator combined with a shift register random integer generator, 
based on the work of Marsaglia.) we obtain for H(l) of 0.47 ±0.01, but also in 
this case H(2) is 0.5 within the errors. This shows that our method is powerful 
and robust and does not suffer of bias as other methods do. On the other hand, 
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the estimations of (3 from the power spectrum have values around 1.8 (instead 
of 2), showing therefore that this other method is affected by a certain bias. 



5 Results and Discussion 

5. 1 Computation of the generalized Hurst exponent 

In this section we report and discuss the results for the scaling exponents H{q) 
computed for q — 1 and q = 2. These exponents H(l) and H(2) for all the as- 
sets and different markets (presented in Section 2) are reported in Figs. 4 and 5 
respectively. Figures 4 (a) and 5 (a) refer to the Treasury and Eurodollar rates 
in the time period from 1990 to 1996. Whereas Figures 4 (b) and 5 (b) are 
relative to the Stock Market indices and Foreign Exchange rates in the time 
period reported in Tables 1 and 2. The data points are the average values of 
H(l) and H(2) computed from a set of values corresponding to different r max 
and the error bars are their standard deviations. Let us first observe that, for 
fixed income instruments (Figs. 4 (a) and 5 (a)), H(2) is close to 0.5 while 
H(l) is rather systematically above 0.5 (with the 3 months Eurodollar rate 
that shows a more pronounced deviation because it is directly influenced by 
the actions of central banks). On the other hand, when Stock markets are 
concerned, we find that the generalized Hurst exponents H(l), H{2) show 
remarkable differences between developed and emerging markets. In particu- 
lar, the values of H(\), plotted in Fig. 4 (b), present a differentiation across 
0.5 with high values of H(l) associated with the emerging markets and low 
values of H(l) associated with developed ones. Moreover, we can see from 
Fig. 5 (b) that the different assets can be classified into three different cat- 
egories: First the ones that have an exponent H{2) > 0.5 which includes all 
indices of the emerging markets and the BCI 30 (Italy), IBEX 35 (Spain) and 
the Hang Seng (Hong Kong). A second category concerns the data exhibiting 
H(2) ~ 0.5 (within the error bars). This category includes: FTSE 100 (UK), 
AEX (Netherlands), DAX (Germany), Swiss Market (Switzerland), Top 30 
Capital (New Zealand), Telaviv 25 (Israel), Seoul Composite (South Korea) 
and Toronto SE 100 (Canada). A third category is associated with H{2) < 0.5 
and includes the following data: Nasdaq 100 (US), S&P500 (US), Nikkei 225 
(Japan), Dow Jones Industrial Average (US), CAC 40 (France) and All Or- 
dinaries (Australia). We find therefore that all the emerging markets have 
H(2) > 0.5 whereas all the well developed have H{2) < 0.5. 

For what concerns the Foreign exchange rates, we find that they show H(l) > 
0.5 quite systematically which is consistent with previous results computed 
with high frequency data [1], although the values here are slightly lower. An 
exception with pronounced H(l) < 0.5 is the HKD/USD (Hong Kong) (Fig. 



9 



4 (b)). This FX rate is or has been at one point pegged to the USD, that is 
why its exponent differs from the others. Whereas in the class H(l) ~ 0.5 
we have: IT L /USD (Italy), PHP/USD (Philippines), AUD/USD (Aus- 
tralia), NZD/USD (New Zealand), ILS/USD (Israel), CAD /USD (Canada), 
SGD/USD (Singapore), NLG/USD (Netherlands) and JPY/USD (Japan). 
On the other hand the values of H(2) (Fig. 5 (b)) show a much larger ten- 
dency to be < 0.5 with some stronger deviations such as: HKD/USD (Hong 
Kong), PHP/USD (Philippines), KRW/USD (South Korea), PEN /USD 
(Peru) and TRL/USD (Turkey). Whereas values of H{2) > 0.5 are found in: 
GBP/USD (United Kingdom), PESO /USD (Mexico), INR/USD (India), 
IDR/USD (Indonesia), TWD/USD (Taiwan) and BRA/USD (Brazil). 



Let us remind that H(2) > 0.5 is commonly associated with a persistent 
behavior in the fluctuations of the returns, whereas an exponent H{2) < 0.5 
indicates anti-persistence [54]. 



These analysis has been done also on different time periods and the values 
are reported in Table 6 for the exponents H(l) and H{2) for the time pe- 
riod from 1997 to 2001 for Foreign Exchange rates, Stock Market indices and 
Treasury rates. In order to verify the stability of these results over different 
time periods we calculate and compare, for the Stock market data, the gen- 
eralized Hurst exponents H(l) and H(2) in the whole time period (shown in 
Table 2) and in time periods of 250 days. Moreover, we tested the numerical 
robustness of our results by using the Jackknife method [60] which consists in 
taking out randomly 1/10 of the sample and iterates the procedure 10 times 
(every time taking out data which were not taken out in the previous runs). 
We observe (see Fig. 6) that the generalized Hurst exponents computed on 
these Jackknife-reduced time series are very close to those computed on the 
entire series with deviations inside the errors estimated by varying T max (as 
described in Section 4). On the other hand, the analysis on sub-periods of 
250 days shows fluctuations that are larger than the previous estimated er- 
rors (and larger than the variations with the Jackknife method) indicating 
therefore that there are physically-significant changes in the market behaviors 
over different time periods (Fig. 6 (a)). This phenomenon was also detected in 
[7] when studying Exchange rates that were part of the European Monetary 
System. It seems that H(l) is particularly sensitive to institutional changes 
on the market. This study confirms it for Stock indices. The scaling exponents 
cannot be assumed to be constant over time if a market is experiencing ma- 
jor institutional changes. Nevertheless, well developed markets have values of 
H(2) that are on average smaller than the emerging ones. Moreover, the weak- 
est markets have oscillation bands that stay above 0.5 whereas the strongest 
have oscillation bands that contain 0.5. 
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5.2 Spectral analysis 



In order to empirically investigate the statistical properties of the time series 
in the frequency domain we perform a spectral analysis computing the power 
spectral density (PSD) [61] by using the periodogram approach, that is cur- 
rently one of the most popular and computationally efficient PSD estimator. 
This is a sensitive way to estimate the limits of the scaling regime of the data 
increments. The results for some SM, FX, TR and ER data in the time pe- 
riods 1997-2001 and 1990-1996, are shown in Figs. 7, 8, 9. For SM and FX 
we compute the power spectra of the logarithm of these time series. As one 
can see the power spectra show clear power law behaviors: S(f) ~ f" 13 . This 
behavior holds for all the other data. 

The non-stationary features have been investigated by varying the window- 
size on which the spectrum is calculated from 100 days to up to the entire 
size of the time series. The power spectra coefficients are calculated through a 
mean square regression in log-log scale. The values reported in Fig. 10 are the 
average of the evaluated f3 over different windows and the error bars are their 
standard deviations. Fig. 10 (a) refers to a time period between 1990 to 1996 
whereas the Stock Market indices and Foreign Exchange rates (Fig. 10 (b)) 
are analyzed over the time periods reported in Tables 1 and 2. Moreover, 
the averaged (3 values in a different time period, namely from 1997 to 2001 
are reported in Table 7 for Foreign Exchange rates, Stock Market indices and 
Treasury rates. These values differ from the spectral density exponent expected 
for a pure Brownian motion (j3 = 2) [62]. However, we have shown in Section 4 
that this method is biased and we have indeed found power spectra exponents 
around 1.8 for random walks using three different random numbers generators. 

It must be noted that, the power spectrum is only a second order statistic 
and its slope is not enough to validate a particular scaling model: it gives only 
partial information about the statistics of the process. 



5.3 Comparison between the generalized Hurst exponent and the power spec- 
tra 

We here compare the behavior of the power spectra S(f) with the function 
pH(2)-i _ according to Eq. 3 - is the scaling behavior expected in the 

frequencies domain for a time series which scales in time with a generalized 
Hurst exponent H(2). We performed such a comparison for all the financial 
data and we report in Figs. 7, 8 those for Foreign Exchange rates and Stock 
Market indices for Thailand and JAPAN (in the time period 1997-2001). In 
Fig. 9 are reported the comparison for the Treasury and the Eurodollar rates 
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having maturity dates 9 = 10 years and 6=1 year respectively. As one can see 
the agreement between the power spectra behavior and the prediction from 
the generalized Hurst analysis is very satisfactory. This result holds also for 
all the other data. Note that the values of 2H(2) + 1 do not in general coincide 
with the values for the power spectral exponents evaluated by means of the 
mean square regression. The method through the generalized Hurst exponent 
appears to be more powerful in catching the scaling behaviour even in the 
frequency domain. 



6 Conclusion 



Scaling behaviors are rather universal across financial markets. By analyzing 
the scaling properties of the g-order moments (Eq. 1) we show that the general- 
ized Hurst exponent H(q) (Eq. 2) is a powerful instrument to characterize and 
differentiate the structure of such a scaling properties. We show that qH(q) has 
a non-linear dependence on q which is a clear signature of deviations from pure 
Brownian motion and other additive models. The empirical analysis across a 
wide variety of stock markets shows that the exponent H(2) is sensitive to the 
degree of development of the market. On one end we find: Nasdaq 100 (US), 
S&P500 (US), Nikkei 225 (Japan), Dow Jones Industrial Average (US), CAC 
40 (France) and All Ordinaries (Australia); all with H(2) < 0.5. Whereas, 
on the opposite side, we find the Russian AK&M, the Indonesian JSXC, the 
Peruvian LSEG, etc. (Fig. 5 (b)); all with H(2) > 0.5. We observe emerging 
structures also in the scaling behaviors of interest rates and exchange rates. 
The robustness of the present empirical approach is tested in several ways: 
by varying the maximum time-step (r max ); by using the Jackknife method; 
by varying the time-window sizes; by comparing with three distinct simulated 
Brownian motions. We verify that the observed differentiation among different 
degrees of market development is clearly emerging well above the numerical 
fluctuations. Finally, from the comparison between the empirical power spec- 
tra and the prediction from the scaling analysis (Eq. 3, Figs. 7, 8 and 9) we 
show that the method through the generalized Hurst exponent describes well 
the scaling behavior even in the frequency domain. 
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Appendix:The Hurst exponent 



Let us consider a time series X(t) defined at discrete time intervals t = 
u, 2u, 3u, ... ku. Let us define the average over a period T (which must 
be an entire multiple of v) as 

Wt = ^!>(H • (5) 
1 k=i 



The difference between the maximum and the minimum values of X(t) in the 
interval [is, T] is called the range R, which is defined as: 

R(T) = max[X(t)] v < t < T - min[X{t)} u < t < T . (6) 



The Hurst exponent H is defined from the scaling property of the ratio: 

R(T ^( T -)" , (7) 



S(T) V v 



where S(T) is the standard deviation: 



S(T) = 



\ 



V 



T/u 



(8) 



k=l 



The Hurst exponent is sensitive to the long-range statistical dependence in the 
signal. It was proved by Hurst [33] and Feller [62] that the asymptotic behavior 
for any independent random process (Poisson process) with finite variance is 
given by: 

* (T > - ^rV /2 (9) 



S(T) \2v 



which implies H = 1/2. However, many processes in nature are not inde- 
pendent random processes, but on the contrary show significant long-term 
correlations. In this case the asymptotic scaling law is modified and R/S is 
asymptotically given by the power law behavior in Eq. 7 with H ^ 0.5. It must 
be noted that the lack of robustness of the original Hurst R/S approach in 
the presence of short memory, heterocedasticity, multiple scale behaviors has 
been largely discussed in the literature (see for instance [34,63,64]) and there- 
fore several alternative approaches have been proposed. Also the fact that the 
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range relies on maxima and minima makes the method error prone because 
any outlier present in the data would have a strong influence on the range. 
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Fig. 1. The Foreign Exchange rates and the Stock Market indices as a function of 
time t in the time period 1997-2001; (a) Nikkei 225; (b) JPY/USD; (c) THB/USD; 
(d) Bangkok SET. 
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Fig. 2. (a) The Treasury rates at 'constant maturity' as a function of t in the time 
period 1997-2001. Each curve corresponds to a maturity date 9, ranging from 3 
months to 30 years and Treasury bill rates to a maturity date 9=3,6 months and 
1 year; (b) The Eurodollar interest rates as a function of t in the time period 
1990-1996. Each curve corresponds to a maturity date 9, ranging from 3 months to 
48 months. 
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Fig. 3. The function qH(q) vs. q in the time period from 1997 to 2001. (a) JAPAN 
(Nikkei 225); (b) JAPAN (JPY/USD); (c) Thailand (Bangkok SET); (d) Thailand 
(THB/USD); (e) Treasury rates having maturity dates 6 = 10 years; (f) Eurodollar 
rates having maturity dates 9 = 1 year. For (f) the time period is 1990 - 1996. 
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Fig. 4. (a) The Hurst exponent H(l) for the Treasury and Eurodollar rates time 
series in the period from 1990 to 1996; (On the x-axis the corresponding maturities 
dates are reported.) (b) The Hurst exponent H (1) for the Stock Market indices and 
Foreign Exchange rates in the time period reported in Tabs. 1 and 2. (On the 
x-axis the corresponding data-sets are reported.) 
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Fig. 5. (a) The Hurst exponent H{2) for the Treasury and Eurodollar rates time 
series in the period from 1990 to 1996; (On the x-axis the corresponding maturities 
dates are reported.) (b) The Hurst exponent H (2) for the Stock Market indices and 
Foreign Exchange rates in the time period reported in Tabs. 1 and 2. (On the 
x-axis the corresponding data-sets are reported.) 
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Fig. 6. (a) The generalized Hurst exponent H(l) for the Stock Market indices in the 
whole time period (see Tab. 2) with its variation (black lines) obtained by using the 
Jackknife method and its variation (dashed lines) when time periods of 250 days are 
considered; (b) The generalized Hurst exponent H (2) for the Stock Market indices 
in the whole time period (see Tab. 2) with its variation (black lines) obtained by 
using the Jackknife method. The square points are the average values of H (1) and 
H(2) computed from a set of values corresponding to different T max . The error bars 
are their standard deviations. 
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Fig. 7. The power spectra of the Foreign Exchange rates compared with the be- 
haviour of y- 2 ^( 2 )~ 1 (straight lines in log-log scale) computed using the Hurst 
exponents values in the time period 1997-2001; (a) Thailand (THB/USD) and (b) 
JAPAN (JPY/USD). The line is the prediction from the generalized Hurst exponent 
H{2) (Eq. 3). 
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Fig. 8. The power spectra of the Stock Market indices compared with the behaviour 
of f- 2H ( 2 )- 1 (straight lines in log-log scale) computed using the Hurst exponents 
values in the time period 1997-2001; (a) Thailand (Bangkok SET) and (b) JAPAN 
(Nikkei 225). The line is the prediction from the generalized Hurst exponent H(2) 
(Eq. 3). 
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Fig. 9. The power spectra compared with the behaviour of f-^i 2 )- 1 (straight lines 
in log-log scale) computed using the Hurst exponents values in the time period 
1997-2001; (a) Treasury rates having maturity dates 6 = 10 years; (b) Eurodollar 
rates having maturity dates 9 = 1 year in the time period 1990-1996. The line is 
the prediction from the generalized Hurst exponent H{2) (Eq. 3). 
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Fig. 10. (a) The averaged (3 values computed from the power spectra (mean square 
regression) of the Treasury and Eurodollar rates time series in the period from 1990 
to 1996; (On the x-axis the corresponding maturities dates are reported.) (b) The 
averaged (3 values computed from the power spectra of the Stock Market indices and 
Foreign Exchange rates in the time period reported in Tabs. 1 and 2. The horizontal 
gray line corresponds to the value of [3 obtained from the simulated random walks 
reported in Table 5. (On the x-axis the corresponding data-sets are reported.) 



Table 1 

Foreign Exchange rates (FX/USD). 



Country 


FX 


Time period 


Hong Kong 


HKD 


1990-2001 


Italy 


ITL 


1993-2001 


Philippines 


rrir 


1001 nnm 

iyyi-zuui 


Australia 


AUD 


1990-2001 


New Zealand 


NZD 


1990-2001 


Israel 


ILS 


1990-2001 


Canada 


CAD 


1993-2001 


Singapore 


SGD 


1990-2001 


Netherlands 


NLG 


1993-2001 


Japan 


JPY 


1990-2001 


Spain 


ESP 


1990-2001 


South Korea 


KRW 


1990-2001 


Hungary 


HUF 


1993-2001 


Germany 


DEM 


1990-2001 


Switzerland 


CHF 


1993-2001 



Country 


FX 


Time period 


United Kingdom 


GBP 


1990-2001 


France 


r rvr 


ryyo-zuui 


Poland 


PLN 


1993-2001 


Peru 


PEN 


1993-2001 


Turkey 


TRL 


1992-2001 


Thailand 


THB 


1990-2001 


Mexico 


PESO 


1993-2001 


Malaysia 


MYR 


1990-2001 


India 


INR 


1990-2001 


Indonesia 


IDR 


1991-2001 


Taiwan 


TWD 


1990-2001 


Russia 


RUB 


1993-2001 


Venezuela 


VEB 


1993-2001 


Brazil 


BRA 


1993-2001 



Table 2 

Stock Market indices (SM). 



Country 


SM 


Time perio< 


United States 


Nasdaq 100 


1990-2001 


United States 


S&P 500 


1987-2001 


Japan 


Nikkei 225 


1990-2001 


United States 


Dow Jones Industrial Average (DJIA) 


1990-2001 


France 


CAC 40 


1993-2001 


Australia 


All Ordinaries (AO) 


1992-2001 


United Kingdom 


FTSE 100 


1990-2001 


Netherlands 


AEX 


1993-2001 


Germany 


DAX 


1990-2001 


Switzerland 


Swiss Market (SM) 


1993-2001 


New Zealand 


Top 30 Capital (T30C) 


1992-2001 


Israel 


Telaviv 25 (T25) 


1992-2001 


South Korea 


Seoul Composite (SC) 


1990-2001 


Canada 


Toronto SE 100 (SE 100) 


1993-2001 


Italy 


BCI 30 


1993-2001 


Spain 


IBEX 35 


1990-2001 


Taiwan 


Taiwan Weighted (TW) 


1990-2001 



Table 2 (continued) 
Country 


SM 


Time 


period 


Argentina 


Merval (ME) 


1993- 


-2001 


Hong Kong 


Hang Seng (HS) 


1990- 


-2001 


India 


Bombay SE Sensex (BSES) 


1990- 


-2001 


Brazil 


Bovespa (BO) 


1993- 


-2001 


Mexico 


Mexico SE (MSE) 


1993- 


-2001 


Singapore 


All Singapore Shared (ASS) 


1990- 


-2001 


Hungary 


Budapest BUX (BUX) 


1993- 


-2001 


Poland 


Wig (WIG) 


1991- 


-2001 


Malaysia 


KLSE Composite (KLSEC) 


1990- 


-2001 


Thailand 


Bangkok SET fBSET") 

U CliAA tlllWIl k_/ J i A \ 1 * V ' 1 J A 1 


1990- 


-2001 


Philippines 


Composite (CO) 


1990- 


-2001 


Venezuela 


Indice de Cap. Bursatil (ICB) 


1993- 


-2001 


Peru 


Lima SE General (LSEG) 


1993- 


-2001 


Indonesia 


JSX Composite (JSXC) 


1990- 


-2001 


Russia 


AK&M Composite (AK&M) 


1993- 


-2001 



Table 3 

Treasury rates (TRi(0)). 
i 9 i 



1 3 months 7 

2 6 months 8 

3 1 year 9 

4 2 years 10 

5 3 years 11 

6 5 years 12 



e 

7 years 
10 years 
30 years 
3 months (Bill) 
6 months (Bill) 
1 year (Bill) 



Table 4 

Eurodollar rates (ERi(9)). 



1 


a 
u 






a 
u 


1 
1 


3 months 


o 


97 


months 


2 


6 months 


10 


30 


months 


3 


9 months 


11 


33 


months 


4 


12 months 


12 


36 


months 


5 


15 months 


13 


39 


months 


6 


18 months 


14 


42 


months 


7 


21 months 


15 


45 


months 


8 


24 months 


16 


48 


months 



Table 5 

Hurst exponents H (1) and H(2) and averaged (3 values computed for random walks 
simulated by using three different random numbers generators: 1) Randn=Normally 
distributed random numbers with mean and variance 1; 2) Rand=Uniformly dis- 
tributed random numbers in the interval (0, 1) and 3) Normrnd=Random numbers 
from the normal distribution with mean and standard deviation 1. These are aver- 
age values on 100 simulations of random walks with 991 and 3118 numbers of data 
points. 

N H(l) H(2) (3 



1) Randn 

991 0.50 ±0.01 0.50 ±0.01 1.8 ±0.1 

3118 0.50 ±0.01 0.50 ±0.01 1.80 ± 0.03 

2) Rand 

991 0.50 ±0.01 0.49 ±0.01 1.8 ±0.1 

3118 0.47 ±0.01 0.50 ±0.01 1.80 ± 0.03 
3) Normrnd 

991 0.49 ±0.01 0.49 ±0.01 1.8 ±0.1 

3118 0.50 ±0.01 0.50 ±0.01 1.80 ± 0.03 



Table 6 

Hurst exponents H(l) and H(2) for Foreign Exchange rates, Stock Market indices 
and Treasury rates in the time period from 1997 to 2001. 



Data 


HQ.) 


H(2) 


Data 


H(l) 


H(2) 


Foreign Exchange rates 












HKD 


0.41 ±0.01 


0.34 ±0.01 


GBP 


0.50 ±0.02 


0.48 ±0.02 


ITL 


0.51 ±0.01 


0.51 ±0.01 


FRF 


0.51 ±0.01 


0.51 ±0.01 


PHP 


0.52 ±0.01 


0.43 ±0.02 


PLN 


0.54 ±0.01 


0.50 ±0.01 


AUD 


0.52 ±0.01 


0.502 ±0.002 


PEN 


0.52 ±0.01 


0.41 ±0.03 


NZD 


0.49 ± 0.01 


0.48 ±0.01 


TRL 


0.56 ±0.01 


0.44 ± 0.04 


ILS 


0.48 ± 0.02 


0.47 ±0.02 


THB 


0.53 ±0.01 


0.50 ±0.02 


CAD 


0.51 ±0.01 


0.48 ±0.01 


PESO 


0.53 ±0.01 


0.50 ±0.01 


SGD 


0.50 ±0.01 


0.47 ±0.03 


MYR 


0.51 ±0.03 


0.45 ± 0.05 


NLG 


0.51 ±0.01 


0.51 ±0.01 


INR 


0.58 ±0.02 


0.53 ±0.01 


JPY 


0.50 ±0.01 


0.49 ±0.01 


IDR 


0.56 ±0.03 


0.53 ±0.03 


ESP 


0.50 ±0.01 


0.49 ±0.01 


TWD 


0.58 ±0.01 


0.51 ±0.01 


KRW 


0.50 ±0.03 


0.39 ±0.06 


RUB 


0.64 ±0.02 


0.47 ±0.03 


HUF 


0.52 ±0.01 


0.52 ±0.01 


VEB 


0.54 ±0.04 


0.49 ± 0.02 


DEM 


0.51 ±0.01 


0.51 ±0.01 


BRA 


0.59 ±0.02 


0.60 ±0.01 


CHF 


0.51 ±0.01 


0.50 ±0.01 









Table 6 (continued) 
Data 


HQ.) 


H{2) 


Data 


H(l) 


H{2) 


Stock Market indices 












Nasdaq 100 


0.47 ±0.01 


0.45 ± 0.01 


TW 


0.53 ±0.01 


0.51 ±0.01 


S&P 500 


0.47 ± 0.02 


0.44 ± 0.01 


ME 


0.57 ±0.01 


0.53 ±0.01 


Nikkei 225 


0.46 ± 0.01 


0.43 ± 0.01 


HS 


0.53 ±0.01 


0.49 ± 0.01 


DJIA 


0.49 ± 0.01 


0.464 ± 0.004 


BSES 


0.54 ±0.01 


0.52 ±0.01 


CAC 40 


0.47 ± 0.02 


0.46 ± 0.02 


BO 


0.51 ±0.01 


0.48 ± 0.01 


AO 


0.49 ± 0.02 


0.46 ± 0.03 


MSE 


0.57 ± 0.01 


0.52 ± 0.01 


FTSE 100 


0.46 ± 0.02 


0.44 ± 0.01 


ASS 


0.57 ±0.01 


0.54 ± 0.02 


AEX 


0.49 ± 0.01 


0.47 ±0.02 


BUX 


0.52 ±0.01 


0.49 ± 0.01 


DAX 


0.50 ±0.01 


0.47 ±0.01 


WIG 


0.49 ± 0.01 


0.44 ± 0.01 


SM 


0.50 ±0.02 


0.48 ± 0.02 


KLSEC 


0.60 ±0.01 


0.51 ±0.02 


T30C 


0.49 ± 0.01 


0.46 ±0.01 


BSET 


0.59 ±0.01 


0.55 ±0.01 


T25 


0.53 ±0.01 


0.51 ±0.01 


CO 


0.59 ±0.01 


0.54 ±0.01 


sc 


0.53 ±0.01 


0.51 ±0.01 


ICB 


0.61 ±0.02 


0.55 ± 0.02 


SE 100 


0.51 ±0.01 


0.48 ±0.01 


LSEG 


0.61 ±0.01 


0.58 ±0.01 


BCI 30 


0.52 ±0.01 


0.48 ± 0.01 


JSXC 


0.57 ±0.02 


0.53 ± 0.02 


IBEX 35 


0.50 ±0.01 


0.48 ± 0.01 


AK&M 


0.65 ± 0.03 


0.51 ±0.01 



Table 6 (continued) 
Data 


HQ.) 


H{2) 


Treasury rates 






TRx 


0.48 ±0.01 


0.44 ± 0.02 


TR 2 


0.55 ±0.01 


0.52 ±0.02 


TR 3 


0.54 ±0.01 


0.52 ± 0.02 


Ti?4 


0.53 ±0.01 


0.52 ±0.02 


TR 5 


0.52 ±0.01 


0.50 ±0.01 


TR 6 


0.51 ±0.02 


0.49 ±0.01 


TR 7 


0.49 ± 0.02 


0.48 ±0.01 


TR 8 


0.52 ±0.01 


0.50 ±0.02 


TR 9 


0.51 ±0.01 


0.48 ±0.01 


TRiq 


0.51 ±0.01 


0.48 ± 0.02 


TR U 


0.56 ±0.01 


0.54 ±0.02 


TR 12 


0.55 ±0.01 


0.53 ± 0.02 



Table 7 

The averaged (3 values computed from the power spectra of the Stock Market indices, 
Foreign Exchange rates and Treasury rates in the time period from 1997 to 2001. 



Data 


Averaged (3 


Data 


Averaged (5 


Foreign Exchange rates 








HKD 


1.6 ±0.2 


GBP 


1.79 ±0.03 


ITL 


1.80 ±0.03 


FRF 


1.81 ±0.04 


PHP 


1.8 ±0.1 


PLN 


1.79 ± 0.04 


AUD 


1.8 ±0.1 


PEN 


1.6 ±0.2 


NZD 


1.8 ±0.1 


TRL 


1.7 ±0.1 


ILS 


1.8 ±0.1 


THB 


1.83 ±0.03 


CAD 


1.80 ±0.03 


PESO 


1.81 ±0.04 


SGD 


1.81 ±0.02 


MYR 


1.8 ±0.1 


NLG 


1.81 ±0.04 


INR 


1.8 ±0.1 


JPY 


1.9 ±0.1 


IDR 


1.83 ±0.04 


ESP 


1.80 ±0.04 


TWD 


1.8 ±0.1 


KRW 


1.8 ±0.1 


RUB 


2.1 ±0.3 


HUF 


1.80 ± 0.03 


VEB 


1.8 ±0.1 


DEM 


1.81 ±0.03 


BRA 


2.0 ±0.2 


CHF 


1.8 ±0.1 







Table 7 (continued) 

Data Averaged (3 Data Averaged (5 



Stock Market indices 








Nasdaq 100 


1.7 ±0.1 


TW 


1.9 ±0.1 


S&P 500 


1.8 ±0.1 


ME 


1.8 ±0.1 


Nikkei 225 


1.8 ±0.1 


HS 


1.8 ±0.1 


DJIA 


1.80 ±0.03 


BSES 


1.82 ±0.03 


CAC 40 


1.8 ±0.1 


BO 


1.80 ±0.02 


AO 


1.8 ±0.1 


MSE 


1.9 ±0.1 


FTSE 100 


1.81 ±0.03 


ASS 


1.9 ±0.1 


AEX 


1.8 ±0.1 


BUX 


1.82 ±0.04 


DAX 


1.8 ± 0.1 


WIG 


1.8 ± 0.1 


SM 


1.8 ± 0.1 


KLSEC 


1.8 ± 0.1 


T30C 


1.8 ±0.1 


BSET 


1.9 ±0.1 


T25 


1.9 ± 0.1 


CO 


2.0 ± 0.2 


SC 


1.9 ± 0.1 


ICB 


2.0 ± 0.2 


SE 100 


1.9 ± 0.1 


LSEG 


2.0 ±0.2 


BCI 30 


1.9 ±0.1 


JSXC 


1.9 ±0.1 


IBEX 35 


1.8 ±0.1 


AK&M 


1.9 ±0.2 


Treasury rates 








TRi 


1.8 ±0.1 


TR 7 


1.9 ±0.1 


TR 2 


1.83 ±0.04 


TR 8 


1.9 ±0.1 


TR 3 


1.86 ±0.05 


TR 9 


1.8 ±0.1 


TR 4 


1.88 ±0.06 


TRiq 


1.82 ±0.04 


TR 5 


1.9 ±0.1 


TR U 


1.85 ±0.04 


TRq 


1.9 ±0.1 


TR 12 


1.9 ±0.1 



